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Abstract — We present a basic high-level structures 
used for developing quantum programming languages. 
The presented structures are commonly used in many 
existing quantum programming languages and we use 
quantum pseudo-code based on QCL quantum program- 
ming language to describe them. 

We also present the implementation of introduced 
structures in GNU Octave language for scientific comput- 
ing. Procedures used in the implementation are available 
as a package quantum-octave, providing a library of 
functions, which facilitates the simulation of quantum 
computing. This package allows also to incorporate high- 
level programming concepts into the simulation in GNU 
Octave and Matlab. As such it connects features unique 
for high-level quantum programming languages, with the 
full palette of efficient computational routines commonly 
available in modern scientific computing systems. 

To present the major features of the described package 
we provide the implementation of selected quantum 
algorithms. We also show how quantum errors can be 
taken into account during the simulation of quantum 
algorithms using quantum-octave package. This is 
possible thanks to the ability to operate on density 
matrices. 

Index Terms — quantum information, quantum pro- 
gramming, models of quantum computation 

I. Introduction 

Quantum information theory aims to harness the 
quantum nature of information carriers in order to 
develop more efficient algorithms and more secure 
communication protocols Q), Q, |0, |4). Unfortu- 
nately, counterintuitive laws of quantum mechanics 
make the development of new quantum information 
processing procedures highly non-trivial task. This can 
be seen as one of the reasons why only few truly 
quantum algorithms were proposed [5|, [6|. 

The laws of quantum mechanics are in many cases 
very different from those we know from the classical 
world. That is why one needs to seek for the novel 
methods for describing information processing which 
involves quantum elements. To this day several formal 
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models were proposed for the description of quantum 
computation process Q, 10, 0, CE), 03], E). 

The most popular of them is the quantum circuit 
model [7], which is tightly connected to the physical 
operations implemented in a laboratory. It allows to 
operate with the basic ingredients of quantum infor- 
mation processing - namely qubits, unitary evolution 
operators and measurements. However, it does not 
provide too much flexibility concerning operations on 
more sophisticated elements required to develop scal- 
able algorithms and protocols eg. quantum registers or 
classical controlling operations. 

Another model widely used in the study of theoret- 
ical aspects of quantum information processing is the 
quantum Turing machine J7J. This model is mainly 
used in the analysis of quantum complexity prob- 
lems [13 1 . Its main advantage it that it provides method 
of comparing efficiency of classical and quantum al- 
gorithms. Unfortunately quantum Turing machine, in 
analogy to its classical counterpart, operates on very 
basic data and thus it cannot be easily used to construct 
quantum algorithms. 

Both quantum circuit model and quantum Turing 
machine share some serious drawbacks concerning 
lack of support for high-level programming and very 
limited flexibility. These problems were addressed in 
the recent research in the area of quantum program- 
ming languages JT4), fl~5), fl6l . 

Quantum programming languages are based on the 
Quantum Random Access Machine (QRAM) model. 
QRAM is equivalent, with respect to its computational 
power, to the quantum circuit model or quantum 
Turing machine. However, it has strictly distinguished 
two parts: quantum and classical. The quantum part 
is responsible for performing parts of a algorithm 
which cannot be computed efficient by a classical 
machine. The classical part is used to control quantum 
computation. This model is used as a basis for most 
quantum programming languages lTl4l . 

Among high-level programming languages designed 
for quantum computers we can distinguish imperative 



and functional languages. The later are seen by many 
researchers as a means of providing robust and scal- 
able methods for developing quantum algorithms 1 17]. 
We, however, focus on the imperative paradigm as it 
provides more straightforward way of implementing 
high-level quantum programming concepts. 

This paper is organized as follows. In Section [III we 
briefly describe the QRAM model of quantum com- 
puter and introduce the quantum pseudocode, which 
was designed to describe this model. In Section [III] 
we introduce high-level programming structures used 
in quantum programming languages. These structure 
are based on the QRAM model of quantum computer. 



In Section IV the implementation of presented con- 
cepts is described and quantum-octave package 
is presented. In Section [V] implementation of quan- 
tum algorithms using quantum-octave package 
is presented. Also the analysis of quantum errors is 
provided in the case of quantum search algorithm. 
Finally Section [Vl] summarize the presented work and 
provides reader with some concluding remarks. 

II. QRAM MODEL OF QUANTUM COMPUTATION 

Quantum random access machine is interesting for 
us since it provides convenient model for developing 
quantum programming languages. However, these lan- 
guages and basic concepts used to develop them are 
our main area of interest. For this reason here we 
provide only the very basic introduction to the QRAM 
model. Detailed description of this model is given in 
|[T8l and |fl9l together with the description of hybrid 
architecture used in quantum programming. 

A. Classical RAM model 

As in many situations in quantum information sci- 
ence, the QRAM models is based on the concepts 
developed to describe classical computational process 
- in this case on the Random Access Machine (RAM) 
model. The classical model of random access machine 
(RAM) is the example of more general register ma- 
chines E<3, ED, El- 

The Random Access Machine consists of an un- 
bounded sequence of memory registers and finite num- 
ber of arithmetic registers. Each register may hold an 
arbitrary integer number. The programme for the RAM 
is a finite sequence of instructions II = (m, . . . , ir n ). 
At each step of execution register i holds an integer 
Ti and the machine executes instruction ir K , where 
k is the value of the programme counter. Arithmetic 
operations are allowed to compute the address of a 
memory register. 

Despite the difference in the construction between 
Turing machine and RAM, it can be easily shown that 
Turing machine can simulate any RAM machine with 
polynomial slow-down only GD . The main advantage 



of the RAM models is its resemblance with the real- 
world computers. 

B. Quantum RAM model and quantum pseudocode 

Quantum Random Access Machine (QRAM) model 
is build as an extension of the classical RAM model. 
Its main goal is to provide the ability to exploit quan- 
tum resources. Moreover, it can be used to perform any 
kind of classical computation. The QRAM allows us 
to control operations performed on quantum registers 
and provides the set of instructions for defining them. 
Schematic presentation of such architecture is provided 
in Figure [T] 
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Fig. 1. The model of classically controlled quantum machine 1 19 |. 
Classical computer is responsible for performing unitary operations 
on quantum memory. The results of quantum computation are 
received in the form of measurement results. 

The quantum part of QRAM model is used to 
generate probability distribution. This is achieved by 
performing measurement on quantum registers. The 
obtained probability distribution has to be analysed 
using a classical computer. 

Quantum algorithms are, in most of the cases, 
described using the mixture of quantum gates, math- 
ematical formulas and classical algorithms. The first 
attempt to provide a uniform method of describing 
quantum algorithms was made in |23l . where the 
author introduces a high-level notation based on the 
notation known from computer science textbooks [24], 



In lfl8l Knill introduced the first formalized lan- 
guage for description of quantum algorithms. More- 
over, it was tightly connected with the model of 
Quantum Random Access Machine. 

Quantum pseudocode proposed by Knill lfT8ll is 
based on conventions for classical pseudocode pro- 
posed in l24l Chapter 1]. Classical pseudocode was 
designed to be readable by professional programmers, 
as well as people who had done a little program- 
ming. Quantum pseudocode introduces operations on 
quantum registers. It also allows to distinguish be- 
tween classical and quantum registers. In quantum 
pseudocode quantum registers are distinguished with 
an underline. They can be introduced by applying 
quantum operations to classical registers or by calling 
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a subroutine which returns a quantum state. In order 
to convert a quantum register into a classical register 
measurement operation has to be performed. 

The example of quantum pseudocode is presented 
in Listing [T] It shows the main advantage of QRAM 
model over quantum circuits model - the ability to 
incorporate classical control into the description of 
quantum algorithm. 
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Fig. 2. Quantum circuit representing quantum Fourier transform 
for three qubits. Elementary gates used in this circuit are described 
in Q]. 



Procedure: F0URiER(a, d) 

Input: A quantum register a with d qubits. Qubits arenumbered 
form to d — 1. 

Output: The amplitudes of a are Fourier transformed over Z 2 d. 
C: assign value to classical variable 

C: perform sequence of gates 
for i — d — 1 to i — 
for j — d — 1 to j = i + 1 
if aj then TZ 2d _i_ 1+ j (a,) 
C: number of loops executing phase depends on 
C: the required accuracy of the procedure 
H(ch) 

C: change the order of qubits 
for j — to j = | - 1 

Listing 1 . Quantum pseudoceode for quantum Fourier transform on 
d qubits. Quantum circuit for this operation with d = 3 is presented 
in Figure [2] 



Operation W(a,i) executes a quantum Hadamard 
gate on a quantum register and SWAT'(ai, dj) 
performs SWAP gate between a,; and cij. Operation 
T?.^ (cti) executes a quantum gate R(4>) is defined as 



R(4>) 



(1) 
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on the quantum register a^. Using conditional construc- 
tion 

if aj then TZ^^ai) 

it is easy to define controlled phase shift gate. Similar 
construction exists in QCL quantum programming lan- 



guage [19|. In Section IV we describe implementation 
of this construction in quantum-octave. 

The measurement of a quantum register can be 
indicated using an assignment 



C. Requirements for quantum programming language 

Taking into account QRAM model we can formulate 
basic requirements which have to be fulfilled by any 
quantum programming language |9|, [26 1, [15|. 

• Completeness: Language must allow to express 
any quantum circuit and thus enable the pro- 
grammer to code every valid quantum programme 
written as a quantum circuit. 

• Extensibility: Language must include, as its sub- 
set, the language implementing some high level 
classical computing paradigm. This is important 
since some parts of quantum algorithms (for ex- 
ample Shor's algorithm) require non-trivial clas- 
sical computation. 

• Separability: Quantum and classical parts of the 
language should be separated. This allows to exe- 
cute any classical computation on purely classical 
machine without using any quantum resources. 

• Expressivity: Language has to provide high level 
elements for facilitating the quantum algorithms 
coding. 

• Independence: The language must be indepen- 
dent from any particular physical implementation 
of a quantum machine. It should be possible to 
compile a given programme for different archi- 
tectures without introducing any changes in its 
source code. 

III. High-level programming structures 
A. Quantum memory 

Quantum memory is a set of qubits indexed by 
integer numbers. Quantum register is a set of in- 
dices pointing to distinct qubits. We will denote those 
registers as r 1 ,r 2 , ■ ■ ■ or in case of single qubits 

as Qi i Q2 j The state of a quantum memory is a 

quantum state of size equal to the number of qubits. In 
the case of quantum-octave we operate on density 
matrices (although some operations on state vectors 
are allowed). We will denote the state of the quantum 
memory by p. 

Following operations on quantum memory are al- 
lowed: 

• Allocation of new register of size n: 



Pt+i =Pt® IP---QXP---QI, 



(2) 



where ® denotes tensor product, | ) the column 

vector and (-| the dual vector. 

Deallocation of a register indexed by register r: 

Pt+i = Tr r (p t ) , (3) 

where Tr r (p) denotes partial trace of p with 
regard to the subsystem indexed by r. 
Unitary evolution U of the quantum memory: 

p t+1 = U Pt U^. (4) 



3 



Application of quantum channel Ki on the quan- 
tum memory: 



pt+i = y^^KiPtKj^. 



(5) 



• Measurement in the computational basis: 

Pt+i=^2\i){i\pt\i){i\, (6) 

i 

P(i)=Tr(\i)(i\p t ), (7) 

where i enumerates the states of computational 
basis. 

For a solid introduction to quantum computation the 
reader may refer to book by Nielsen and Chuang (TJ, 
where all the needed notions are explained in detail. 

In quantum computation, construction of the uni- 
tary gate is the essential part of quantum algorithm 
(program) design process. It is a difficult task to 
write a quantum program using only elementary set 
of gates ie. CNot and one qubit rotations. Therefore it 
is desirable to introduce some techniques that facilitate 
the process of composition of complex quantum gates. 
Some of those techniques are presented below. We 
will refer to implementation of those techniques in 
quantum-octave which is described in details in 
section llVl 

B. Composed and controlled gates 

1) Composed gate: Given one-qubit unitary gate 
G, quantum register r, and size of the gate s we can 
construct composed gate U" according to the formula: 



i Xi , where X; 



G if i e r, 
I if i i r 



(8) 



2) Controlled gate with multiple controls: Given 
one-qubit unitary gate G, quantum register r c we call 
control, and quantum register r t we call target, and 
size of the gate s we can construct controlled gate 
U", r according to the formula: 

U S n \r c = ®Ll X i + ®Ll Y h Where 

' |0)(0| if i£r c , 

I if i <£ r c , 

G if i e r t , 

Y: ^ : { |1)(1| if ier c , 

I if i 4 r c U r t 



Xi 



■ (9) 



We assume that r t n r c = 0. Sometimes we will omit 
the size parameter s. 

C. Conditionals 

One of high-level technique used in quantum pro- 
gramming are quantum conditions [27 1. The main idea 
behind quantum conditions is construction of quantum 
gates controlled by predicates on control registers. 



1 ) Condition on quantum variable: The if-then-else 
structure controlled by a quantum variable and acting 
on a quantum variable was introduced in QCL [28 1, 

na. 

In Figure [3] the realisation of this concept is pre- 
sented. If qubit go is in the state |1) the G\ gate is 
applied to qubit qi, otherwise the gate Gi is applied. 

We may write this circuit in the following way: 

IF qo (G iqi )ELSE(G 2qi ) = Not qo G 2qilqo Not qo G lqilqo . 

(10) 



qbit gi,ij2 

if (qi) then 

Si (<? 2 ) 
else 

62(92 ) 



ki) - Gi 



— 
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Fig. 3. Example of a simple quantum if-then-else structure. 

For a given control register r c , target register r t 
and two quantum gates G± and G2, we may define 
quantum condition in the more general way, namely 

IF rc (G irt )ELSE(G 2 r t )= (11) 
Yl (NotiG 2rt \ rc Noti) G Xrt \r c , 

ieV(r t )\{0} 

where V(-) denotes the power set. 

2) Condition on mixture of classical and quantum 
variables: One may consider the relation between the 
state of quantum register and the value of the classical 
variable. By [[x]] r we will denote a numerical value of 
ordered (in ascending order) elements of the register 
x in regard to register r, for example the value of 
[[{4, 9}]]{ 2 ,4,7,9} is 10. By [r] we will denote the 
value of the register in order to use it as argument for 
arithmetic comparison. For example [r] < 4 means: 
"all those values of r that are less than four." 

Code and circuit in Figure |4] show the idea and the 
implementation of conditional operation controlled by 
expression 'less than' operating on classical constant 
and quantum register. 



Less than 



pseudocode 



qnibble r 
qbit 91,92 

if (r < 4) then 

61(92) 
else 

62(91) 



quantum-octave 



r = newregister (4 
ql=newregister ( 
q2=newregister ( 



qif (... 

qrlt(qureg(ql) 
{GI , qureg (q2) ) 
{G2, qureg (ql ) ) 

Quantum Circuit 



); 
i) ; 
i) ; 



.4) 




Fig. 4. Example of quantum conditional operation with inequality. 
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In the general case, gate implementing any relation 
(marked as ©) can be constructed in the follow- 
ing way: 

IF [rc]mN {G 1 r t )ELSE{G 2 rt) = (12) 



[r c }®N 

Yl (NoUG 2rt \ r NoU) J] (NotiG 
ieT ieT 
where sets T and T are defined as follows 



lr t \r c Noti) , 



T = T(r c ) \ {x\x e T(r c ) A [[x]] P „ © N}, (13) 
T = V(r c ) \ {x\x e T(r c ) A -.([[x]] r „ © N)}. (14) 

Note that T U T = V{r c ). 

In quantum-octave standard arithmetic relations 
=,7^, <,>,<,> are implemented. 

D. Expressions 

We may consider more complicated expression on 
quantum registers. For example logical operators and 
quantum pointers. Logical operators allow to apply 
an controlled operation to the target register only if 
a given logical expression on control registers is true. 
A quantum pointer allows to apply controlled gate on 
the target register selected by the state of the control 
register. 

1) Logical expressions on quantum variables: The 
gate that implements logical expression (denoted here 
by o) is constructed according to the following equa- 
tion: 

IF [r cl ]®i Wio[r C2 ]©2JV2 (G 1 r t )ELSE(G 2 r t ) = 
= l[ (NattG 2 r t \r B Noti) JJ (NotiG ln \ re Noti) , 



ieT 



ieT 



(15) 



where sets T and T are defined as follows: 

T = V(r c ) \ {xi U x 2 \xx C r Cl ,x 2 C r C2 (16) 

A ([[Sl]]r ei ©1 Ni O [[X 2 }] rc2 © 2 A^ 2 )}, 

T = V(r c ) \ {xi U x 2 \xx C r Cll x 2 C r C2 (17) 
A- ([[aJi]]r ai ©i Ni o [[X2]] ra3 ©2 ^2)} 

and r c = r Cl U r C2 . 

An example of quantum conditional gate controlled 
by logical expression defined on quantum registers is 
presented in Figures [5] and [6] 

2) Quantum pointers: In analogy to the concept of 
pointers and indirect addressing in classical program- 
ming, one may introduce quantum pointers. The idea 
is to use control register to control on which of the 
target registers an operation should be applied. 

Let us assume that one has the n-bit control register 
and a set of 2" -bit target registers. The control register 
stores the address of target register to which given 
unitary operation shall be applied. In order to visualise 
the use of a quantum pointer an example is shown in 
Figure [7] 



And 



pseudocode 



quantum- octave 



qbit 91,92,93,94 



if ( 
Si 
else 



91 and 

(-33) 

(-34) 



92) then 



ql = newregister (1) 
q2=newregister (1) 
q3=newregister(l) 
q4= ne wr e gi s ter ( 1 ) 

qif (... 
qrand ( . . . 
qreq(qureg(ql) , 
qreq(qureg(q2) , 
{Gl , qureg (q3) ) , 
{G2,qureg(q4) )) 



Quantum Circuit 
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Fig. 5. Example of quantum conditional operation with "and" 
operator. 



Or 



pseudocode 



qbit 91,92,93,94 

if (91 or 92) then 

61(93) 
else 

62(94) 



quantum- octave 

ql = newregister(l) 
q2=newregister ( 1 ) 
q3 = newregister(l) 
q4=newregister(l) 

qif (... 
qror ( . . . 

qreq(qureg(ql) ,1) ,.. 

qreq(qureg(q2) ,1) ) 
(Gl , qureg(q3) ) , ... 
|G2,qureg(q4) }) 



Quantum Circuit 
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Fig. 6. Example of quantum conditional operation with "or" 
operator. 



Formally, a quantum pointer controlled by register 
r c with target r t is constructed in the following way: 



POINT rt {G [rc] )= H (Not rc \iG 



Mr Jr. 



Not 



ieV(r c ) 



r a \i) 



(18) 



IV. Package quantum-octave 

Package quantum-octave ll29l . POl . Il3ll pro- 
vides a quantum programming, simulation and analysis 
language implemented as a library of functions for 
GNU Octave 11321 . 

GNU Octave is computer algebra system (CAS) and 
high level programming language designed primarily 
to perform numerical calculations. The basic data 



5 



Pointer 



pseudocode 


quantum-octave 


qreg[2] gi 
qnibble q-z 

if (* qi ) then 

Q(Q2) 


q 1 = n e w register (2) ; 
q2=newregister(4) ; 

qpointer (G, qureg (ql ) ,qureg(q2) ) 



Quantum Circuit 
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Fig. 7. Example of simple quantum conditional operation controlled 
by quantum pointer. 

structure in Octave is a matrix (integer, real or com- 
plex), therefore it is a natural choice as a basis for the 
implementation of quantum programming language. 

GNU Octave supports sparse matrices and dis- 
tributed computing in shared and distributed memory 
models. GNU Octave is a very flexible and easily 
extendible tool. It is also free software and it can be 
used in a wide range of operating systems. 

A. Design choices 

The main goal of the design of quantum-octave 
is to provide a flexible and useful tool for simulation of 
quantum information processing. Therefore it is based 
on GNU Octave a high-level scientific programming 
language. This allows for a seamless integration of 
very efficient matrix operations and numerical proce- 
dures with the library of specialized functions provided 
by quantum-octave. As GNU Octave is to large 
extent compatible with Matlab, provided functions 
can be also used to simulate and analyse quantum 
algorithms in Matlab. 

One of the unique features of quantum-octave 
is its ability to operate on both pure and mixed quan- 
tum states. It allows to perform unitary as well as non- 
unitary evolution represented by quantum channels. 

Quantum gates can be constructed by the user in 
various ways: by calling provided subroutines, by 
building their own subroutines, by using quantum 
control structures. Additionally the user can build and 
use quantum channels or use those already provided. 
Most of the quantum-octave functions operate on 
quantum registers and therefore the quantum opera- 
tions build with their use are re-allocable. 

A good illustration of those features is presented in 
the following listing [2]that contains the implementation 



of Quantum Fourier transform in quantum-octave. 
It can be compared to pseudo code version of the same 
procedure listed in Listing [T] 

function ret = dft ( gatesize ) 
2 n= gatesize; 

cir=id (n) ; 
4 for i =[ 1 : n] 

for j=[l:i-l] 
6 cir = circuit(cir,cphase(pi/(2 A (i— j)),j,i, 

gatesize)); 

endfor 

8 cir=circuit(cir,productgate(h,i, gatesize)); 

endfor 

10 ret=cir = circuit(cir , flip(n)); 
endfunction 

Listing 2. Quantum Fourier transform in quantum-octave 

Package quantum-octave can operate on sparse 
and full matrices depending on users choice. Sparse 
matrices need much less memory to store but opera- 
tions on them may be slower. Full matrices tend to 
consume huge memory space, but operations on them 
are generally faster. In case of full matrices it should be 
possible to operate on states of size up to ten qubits on 
a contemporary workstation. Sparse matrices should 
allow to simulate the quantum systems of up to 20 
qubits. 

Although quantum-octave is not, strictly speak- 
ing, a programming language ready to program real 
quantum devices, with some effort it can be trans- 
formed in such a way that it would be able to compile 
a high level programs to some sort of quantum assem- 
bler. One should note that broad range of functions 
allowing the analysis of states is implemented in the 
package. Those function are described in the following 
section. 

B. Description 

Package quantum-octave is designed to allow 
the user to operate on different levels of abstraction. 
User can prepare complex gates and quantum channels 
from basic primitives such as single qubit rotations, 
controlled gates, single qubit channels. Most of the 
functions that form this library operate on quantum 
registers, which makes the preparation of quantum 
gates and channels very "natural". The library is 
implemented in such a way that depending of user's 
choice it may operate on full or sparse matrices. 

quantum-octave can work in two modes: as a 
library or as a programming language/simulator. Li- 
brary mode is the default. To switch to language/simu- 
lator mode one has call quantum_octave_init ( ) 
function. In case of the second mode quan- 
tum-octave allocates and manages an internal 
quantum state and maintains the quantum registers. 
Such functions as evolve (), applychannel ( ) , 
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measurecompbasis () operate directly on the in- 
ternal state. Listings of Deustch's [3] and Grover's [4] 
algorithms show use of language/simulation mode. 

1) Convention: Following conventions are used in 
quantum-octave. 

- Quantum register is horizontal integer vector con- 
taining indices of qubits starting from one. 

- Ket is vertical complex vector. 

- Bra is horizontal complex vector. 

- Density matrix is complex square matrix always 
of dimension n-th power of two by n-th power 
of two. 

- Binary string is 0, 1 -horizontal vector, that en- 
codes a binary number. Order of bits is from MSB 
to LSB. 

- Size of the gate or channel is always given in 
terms of a number of qubits it acts on. If size is 
written in square brackets it means that it can be 
omitted if the gate or channel acts on the whole 
system and quantum-octave was initialised. 

2) Quantum gates: Package quantum-octave 
supplies set of elementary gates known in quantum 
computation. 

- sx, sy, sz - return one-qubit Pauli operators sx 

- o~ x , sy — (J y , sz — er 2 . 

- id(n) - returns identity matrix: I„. 

- roty(a), rotz(a), rotx(a) - return rota- 
tion matrix by angle a around appropriate axis. 

- qft (n) - returns quantum Fourier transform on 
n qubits. 

- swap (size, qubits) - returns swap gate of 
a given size that swaps qubits given as two- 
element vector. 

- qubitpermutation (permutation) 

- returns unitary gate that performs given 
permutation. 

- h - returns one-qubit Hadamard gate. 

- phase (pO , pi ) - returns one-qubit phase gate, 
with p0,pl phase parameters. 

3) Basic functions: Following functions are essen- 
tial to prepare a quantum state and to implement 
a quantum algorithm, protocol or game. 

- ket (binvec) - returns ket for given binary 
string. 

- ketn (int, size) - returns a ket of size 2 slze 
for given integer number. 

- state (pure_state ) - returns density matrix 
for a given ket. 

- mixstates (al , mixed_statel , [a2, 
mixed_state2 , . . . ] ) - returns convex 
combination of density matrices with coefficients 
al , a2 , .... 

- product gate (gate, targetreg [ , size] ) 

- returns a controlled gate of a given size that 
applies given gate on target register. See 



Eq.0 

- controlledgate (gate, controlreg, 
targetreg [, size ] ) - returns a controlled 
gate of given size that applies gate on spec- 
ified target register and is controlled by 
control register. See Eq. [9] 

4) Quantum conditional operations: The functions 
listed below implement quantum conditional opera- 
tions, quantum expressions and pointers. They are 
useful to simplify the implementation. 

- qif (expression, ifpart, elsepart, 
size) - returns quantum gate of given 
size, controlled by expression that ap- 
plies ifpart if expression is true and 
elsepart if expression is false, ifpart and 
elsepart are cellarrays in the form: {gate, 
target_register } . See Eq. [TTj 

- qreq (register, integer ) - returns ex- 
pression: [register] equals integer. See 



Eq. 12 



qrne (register, integer) - returns ex- 
pression: [register] not equal integer, 
qrge (register, integer) - returns ex- 
pression: [register] is greater or equal to 
integer. 

qrgt (register, integer) - returns ex- 
pression: [register] is greater than integer, 
qrle (register, integer) - returns ex- 
pression: [register] is lesser or equal to 
integer. 

qrlt (register, integer) - returns ex- 
pression: [register] is lesser than integer, 
qrin (register, set) - returns expression: 
[register] is in set. 

qror (exprl, expr2) - returns logical or on 



expressions exprl and expr2. See Eq. 15 

- qrand (exprl , expr2 ) - returns logical and 
on expressions exprl and expr2. 

- qpointer (gate, contrregister, 
targteregister [ , size] ) - returns 
quantum gate of given size, controlled by 
controll register that applies gate on 
target register. See Eq. [T8| 

5) Evolution, channels and measurement: The fol- 
lowing group of functions allows to control the evolu- 
tion of quantum states and introduces the application 
of channels and measurement. 

- evolve (evolution [, state ] ) - applies 
unitary evolution to the state, returns the 
result of the evolution. See Eq. [4] 

- channel (name,p) - returns Kraus operators 
acting on one qubit, parametrised by p 
allowed names are: "depolarizing", 
"amplitudedamping", "phasedamping", 
"bitflip", "phaseflip" and 
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"bit phase flip". 

- localchannel (kraus, targetreg [, 
chsize] ) - returns a channel being the 
extension of defined by kraus operators 
channel, acting on target register. 

- applychannel (elements [, state] ) - 
applies on the state non unitary evolution 
defined by set of Kraus operators (elements), 
returns the result of the evolution. See Eq. [5] 

- ptrace (state, targetreg) - returns re- 
duced density matrix for the state with 
target register traced out. See Eq. [3] 

- circuit (gate [ , gate]) - returns circuit 
composed of the sequence of gates. 

- measurecompbasis ([ state ] ) - returns 
the probability distribution of the o z measurement 
on the given state. 

- isunitary (gate) - returns true if the 
gate is unitary, otherwise returns false. 

- ischannel (operators ) - returns true if 
the operators form valid quantum channel, 
otherwise returns false. 

- collapse (distribution) - chooses and 
returns a basis state at random according to 
distribution. 

6) Computation and control: Following functions 
allows to control the quantum heap and configure the 
behaviour of the library. 

- quantum_octave_init ( ) - initialises the 
simulated system, creates quantum state with zero 
qubits allocated and empty list of registers. 

- set_quantum_octave_sparse ( [true | 
false] ) - switches on or off use of sparse 
matrices by all quantum-octave functions. 

- newregister (size) - creates new register 
of given size, allocates qubits on quantum heap, 
returns register id. 

- clearregister (regid) - removes regid 
register from quantum heap. Traces out appropri- 
ate qubits from the internal state. 

- qureg (regid) - returns quantum register 
to which regid points. 

- getstateO - returns the internal quantum 
state. 

7) Well known states: Some of the states commonly 
used in quantum algorithms are implemented in the 
library as separate functions. 

- ghz (n) - returns Greenberger-Horne-Zeilinger 
state for n qubits: ^(|0)® n + |l))® n . 

- phip -returns Bell |$+) state: 4=(|00) + 

- phim - returns Bell |$ _ ) state: ^g(|00) - 

- p sip -returns Bell |#+) state: ^=(|01) + 110}). 

- psim - returns Bell |^~) state: ^(|01) - |10)). 

- maximallymixed (n) - return density matrix 
maximally mixed state: -I n . 



- wernersinglet (a) - returns 2-qubit Werner 
state: a(|00) - |11>)((00| - (11|) + (1 - a)f . 

8) Analysis: Package quantum-octave pro- 
vides standard functions for analysis of quantum states, 
widely used in quantum information literature. Among 
them the most important are: 

- negativity ( state, qubits) - computes 
negativity of the state in respect to qubits. 

- entropy (state) - computes Von Neuman 
entropy of the state. 

- concurrence ( state ) - computes concur- 
rence of the state. 

- fidelity (rho, sigma) - computes fidelity 
between density matrices rho and sigma. 

- f idelitypuremixed (psi, rho) - com- 
putes fidelity between ket psi and density matrix 
sigma. 

- tracenorm ( state) - computes trace norm of 
the state. 

- partialtranspose (state, 
targetreg) - returns matrix being partial 
transposition of state matrix in regard 
to target register. 

The next section presents the applications of quan- 
tum-octave and various programming techniques 
for solutions of quantum programming problems. 

V. Examples and applications 

In what follows the applications of quan- 
tum-octave and various high-level programming 
techniques are discussed . It is shown how quantum 
processes, such as algorithms may be implemented, 
simulated and analysed with this tool. 

A. Deutsch's problem 

One of the simplest quantum algorithms is 
Deutsch's algorithm. Although it may seem trivial, 
this algorithm shows two very important features of 
quantum computation. The first one is taht by taking 
the advantage of a superposition one can compute any 
binary function for all its arguments in one step. The 
second features is is that it is only possible to retrieve 
information about property of a function and not on 
its values. 

Let us assume that we have a black box that is 
usually called the oracle. This box computes a function 
/ : {0, 1} -> {0, 1}. We do not know if that function 
is constant /(0) = /(l) or injective /(0) = /(l). In 
classical case we have to ask the oracle twice to check 
which kind the function / is. But in quantum case it is 
possible to solve this problem asking the oracle only 
once. 

The algorithm goes as follows: 

1) Prepare the state: |*) = |0) <£> |1),. 



8 



2) Apply the Hadamard H® 2 gate on the state |<J>) 
you will get 

|o> + |i>,jo> 



*i = 



i 



V2 



(19) 



V2 ■ 

3) Apply the gate U f : \x)®\y) -> \x)®\f(x)®y) 
on the state you will get: 

' -M±lli M) for constant /, 

for injective /. 



1*2 



± 



%/5 

± \Q)-K 

(20) 

on the state \ you will get: 




for constant /, 

for injective /. 

(21) 

5) Measure state of the first qubit, you will get 
|0) in case of constant function, |1) for injective 
function. 

Quantum circuit representation of Deutsch's algo- 
rithm is presented in Figure [8] The U f gate provide 
a reversible implementation of function / and the 
symbol denotes a measurement. 



|0> ~H 



X 



H 



Uf 



x 



H 



y © f{x) 



Fig. 8. Deutsch's algorithm 

The implementation of Deutsch's algorithm pre- 
sented in listing [3] is an introductory example of 
application of quantum-octave for simulation of 
a quantum algorithm with all basic steps of compu- 
tation: initialization of the quantum computer, unitary 
evolution and measurement. 

Below we describe simulation steps (compare with 
circuit in Figure [8]). Numbers on the left refer to the 
line numbers in Listing [3] 

5 : initialisation of the simulator, 

7, 8 : allocation of registers, 

10-17 : definition of all four possible oracles, 

19 : application of Not on second qubit, 

20 : application of H ® H, 

21 : application of the oracle, 

22 : application of H (g) I, 

24 : tracing out of second register, 
26 : return the probability distribution of the 
measurement outcome. 

l # input : identifier of the function 
# output: state after execution of Deutsch's 
algorithm 
3 function ret = deutsch (num) 

# initialize the simulation 
5 quantum_octave_init () ; 



# declare and allocate registers 
7 rl = newregister ( 1 ) ; 

r2= newr e g i s ter ( 1 ) ; 
9 # declare functions 

f (l}=id(2) ; 
n f {2)= product gate ( sx , qureg ( r2) ) ; 

f ( 3 ( = q i f (qreq(qureg(rl ) ,1) 
13 { sx , qureg ( i"2 ) } , . . . 

( id , qureg(r2) )) ; 
15 f { 4} = qif ( qreq ( qureg ( r 1 ) ,0) 

{ sx , qureg ( r2 ) ) , . . . 
n { id , qureg (r2 ) )) ; 

# do the algorithm 

19 evolve ( productgate ( sx , qureg (r2 ) ) ) ; 

evolve( productgate (h , [ qureg(rl ) ,qureg(r2) ]) ) ; 
21 evolve(f {num } ) ; 

evolve ( productgate (h , qureg ( r 1 ) ) ) ; 
23 # throw away second register 

clearregister (qureg(r2)) ; 
25 # return the outcome 

ret=measurecompbasis () ; 
27 cndfunction 

Listing 3. Deutsch algorithm in quantum-octave 



B. Grover's algorithm 

To illustrate more advanced usage of the presented 
concepts we use the quantum algorithm for searching 
a unordered database. The algorithms was proposed by 
Grover 031 . 04ll . If35l and its detailed description and 
analysis can be found in [36], [37|. Here we present the 
implementation of Grover's algorithm which presents 
the features of quantum-octave related to the ob- 
servation of quantum errors. We show the propagation 
of initial errors during the execution of the algorithm. 

Grover's search algorithm is one of the most im- 
portant quantum algorithms. This is especially true 
since many algorithmic problems can be reduced to 
exhaustive search. However, like in the case of any 
quantum procedure, the efficiency of the algorithm 
depends on the ability to avoid errors during the 
procedure. Thus, it is important how quantum errors 
affect the executions of the algorithm. 

1 ) Statement of the problem: Let X be a set and 
let / : X -> {0, 1}, such that 



m = 



1 x = xq 
43- x 7^ xq 



,x e x, 



(22) 



for some marked xq € X. 

For the simplicity we assume that X is a set of 
binary strings of length n. Therefore \X\ = 2 n and 
/:{0,1}"->{0,1}. 

We can map the set X to the set of states over %® n 
in the natural way as 



x «-> \x). 



(23) 



The goal of the algorithm is to find the marked 
element. This is achieved by the amplification of the 
appropriate amplitude ll36l . Il37l . 

2) The algorithm: The Grover's algorithm is com- 
posed of two main procedures: the oracle and the 
diffusion. 
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a) Oracle: By oracle we understand a function 
that marks one defined element. In the case of this 
algorithm, the marking of the element is done by 
negation of the amplitude of the state that we search 
for. 

With the use of elementary quantum gates the oracle 
can be constructed using an ancilla \q) in the following 

way: 

0\x)\q) = \x)\q®f(x)). (24) 



If the register \q) is prepared in the state: 

|o)-|i) 



|9) = -ff|l) 



V2 



(25) 



then by substitution, equation [24] is re-transformed to: 

-ID . , 10) -II) 



0\x) 



(-l) f W\x) 



and by tracing out the ancilla we get: 

0\x) =-(-l) /(x) |x). 



(26) 



(27) 



Thus the oracle marks a given state by inverting its 
amplitude. 

b) Diffusion: The operator D rotates any state 
around the state 

m = 4^5>>. (28) 

V Z x=0 

D may written in the following form: 

D = -H® n (2\0){0\-I)H® n = 2\i/)){tp\ - 1. (29) 

c) Grover iteration: The first step of the algo- 
rithm is to apply Hadamard gate H® n on all the qubits. 
Then we apply gate G = DO several times. 

d) Number of iterations: Application of diffusion 
operator on the base state \n) gives 

2 



- H® n hH vn \n) = -\x ) 



N 



E 



(30) 



Application of this operator on any state gives 
D\x) = Y^u^x) + ~yY,\v)) 



^2(-a x +2s)\x), 

X 

s = ^E«* 



where 



(31) 
2 n_ 



is arithmetic mean of coefficients a x , x = 0, 
1. 

fc-fold application of Graver's iteration G on initial 
state Is) leads to [361: 



(32) 



G k \s) = a k \x)+Pk\xo), 

x=£x 



with real coefficients 
1 



cos (2k + 1) 9, I3 k = sin (2k + 1) i 



where 9 is an angle that fulfils the relation: 

1 

w 



(33) 



(34) 



Therefore the coefficients ctk,Pk are periodic functions 
of k. After several iterations amplitude of /3*. rises and 
others drop. The influence of the marked state |aro) 
on the state of the register is that the initial state \s) 
evolves towards the marked state. 

The /3k attains its maximum after approximately 
jvN steps. Then it begins to fall. Thus, the number 
of steps needed to transfer the initial state towards the 
marked state is of 0(y/n). In the classical case the 
number of steps is of 0(n). 

e) Measurement: The last step of the Graver's 
algorithm is the measurement. Probability of obtaining 
of the proper result is |/3/c| 2 . 
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Fig. 9. The circuit for Graver's algorithm. 

3) Graphical interpretation: There a exists very 
nice graphical interpretation of Graver's algorithm. 

Let | a) denotes the sum of states orthogonal to the 
state we are searching for |cco) 



1 



=, E i*>. 



(35) 



x^x 



and for consistence we will write \/3) — \xq). Then, 
on the plane spanned by \a) and \(3), we can observe 
of evolution of the state vector. 

By putting values from equation [33] into equation 
32 we get following relation: 



G k \s) =cos((2* + l)0)|a)+ffln((2A + 1)0)|0). (36) 
Exemplar behaviour of this equation for 2 3 states is 



presented in Figure 10 



4) Implementation: Listing [4] presents the imple- 
mentation of function grover. We will apply quan- 
tum noise at the end of each Grover iteration and 
observe its influence on its efficiency. 

To simulate this behaviour we will insert the code 
from Listing [5] after line 21 of the implementation. 
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Fig. 10. Visualisation of Grover's algorithm Projection on 
plane spanned by \a) and Vector \ip) is flat superposition of 
all the possible states. 



i # function implementing Graver ' s algorithm 

# input : number we are looking for , size of the 

system 

3 # output : probability distribution after 

# execution of the algorithm 
5 function ret = gr o v er ( num , S ) 

# initialize the simulation 
7 quantum_octave_init(); 

# allocate register 
9 rl = newregister(S); 

# number of elements 
n N=2 A length(qureg(rl )) 

# calculate number of iterations 
13 k = floor ((pi/4)*sqrt(N)) ; 

# prepare the system in flat superposition of base 

states 

15 evolve ( product gate (h , qureg ( rl ))) ; 

# Grove r iterations 
n for i = l:k 

# ask the oracle 

19 evolve(oracle (num , qureg(rl))); 

# diffu s e 

21 e vol ve ( diffuse (qureg(rl))); 
endfor 

23 # return probability distribution of base states 

ret = measurecompbasis() ; 
25 endfunction 



27 # function implementing oracle 

# input : number to mark , size of the system 
29# output : gate implementing oracle of the size 

function ret = oracle ( num , register) 
31 l = length (register); 

ret = id ( 1 ) ; 
33 ret (num+ 1 ,num+ 1 ) = — 1 ; 

endfunction 

35 

# function implementing oracle 

37 # input : register on which implement diffu sion 

# output : diffu sion gate of the size 2 A I 
39 function ret = diffuse(register) 

l=length ( register ) ; 
41 ret = circuit(... 

pr oductgate ( h , r e g i s te r ,1) 



2 A / 



» (2*ketn(0,l)*bran(0,l) - id ( 1) ) 

p r oductgate (h, register , 1) . . . 

i ); 

endfunction 

Listing 4. Grover's algorithm in quantum-octave 

apply channel ( 
i localchannel( 

channel (channel name ,p) , qureg (rl ) 

i ) 

); 

Listing 5. Adding noise to Grover's algorithm 

a) Simulation results: The results of the simula- 
tion of noisy Grover's algorithm acting on system of 
size from three to six qubits when system is affected by 
noise modelled with depolarizing channel are shown 



in Figure 1 1 One may observe that rate of successful 
application of the algorithm drops quickly with raising 
amount of noise. This effect is more significant for 
larger systems. This result clearly indicates that it 
is not possible to successfully implement Grover's 
algorithm in presence of large amounts of noise if no 
error correction scheme is applied. 



depolarizing channel 




noise amount a 



Fig. 11. Influence of depolarizing channel parametrized by single 
real number a on probability of successful finding of sought element 
in Grover's algorithm implemented with 3, 4, 5 and 6 qubits. 



VI. Summary 

We have introduced an original solution to the 
problem of simulation of quantum processes. This 
solution is provided by quantum-octave, a library 
that is build upon GNU Octave high level program- 
ming language, which provides high-level quantum 
programming structures. 

Although strictly speaking, quantum-octave is 
not a programming language but a library, together 
with GNU Octave, it is very convenient and flexible 
tool. Programs written in quantum programming lan- 
guages, such as QCL, can be easily rewritten using 
this library, thanks to the use of quantum memory, 
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registers and routines. Scalable programs can be eas- 
ily implemented in quantum-octave, so the pro- 
grammer does not have to think about details of the 
implementation. 
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